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Interactions between light and conducting nanostructures can result in a variety of novel and fascinating 
phenomena. These properties may have wide applications, but their underlying mechanisms have not been 
completely understood. From calculations of surface charge density waves on conducting gratings and by 
comparing them with classical surface plasmons, we revealed a general yet concrete picture about coupling of 
light to free electron oscillation on structured conducting surfaces that can lead to oscillating subwavelength 
charge patterns (i.e., spoof surface plasmons but without the dispersion property of classical surface plasmons). 
New wavelets emitted from these light sources then destructively interfere to form evanescent waves. This 
principle, usually combined with other mechanisms (e.g. resonance), is mainly a geometrical effect that can be 
universally involved in light scattering from all periodic and nonperiodic structures containing free electrons, 
including perfect conductors. The spoof surface plasmon picture may provide clear guidelines for developing 
metamaterial-based nano-optical devices. 

PACS numbers: 42.25.Bs, 42.79.Dj, 73.20.Mf, 84.40.Az 



I. INTRODUCTION 

The various novel and unusual optical properties of con- 
ducting nanostructures, such as anomalous diffraction from 
metallic gratings, enhanced light transmission through sub- 
wavelength slits or holes, light polarizing through wire grid 
polarizers, surface-enhanced Raman scattering, negative re- 
fraction imaging, etc, have attracted tremendous attention in 
recent years . 1 ^ 3,4 - 5 To date the coupling of light with surface 
plasmons (SPs) have been widely adopted to explain these 
anomalous phenomena. However, the SP picture elaborated 
in numerous case studies in the literature actually corresponds 
to a very general concept about coupling of electromagnetic 
(EM) waves to free electron oscillation on conducting sur- 
faces that can generate evanescent EM wave modes. This big 
picture is correct without doubt, but it is too general for one 
to obtain a clear and straightforward understanding of the es- 
sential underlying mechanism. Due to this uncertainty, the 
SP-like wave modes have been usually assumed to be sim- 
ilar to the classical SPs (CSPs) on planar metal surfaces, 6 
but this assumption is obviously challenged by the fact that 
(nearly) perfectly conducting structures that do not support 
CSPs still have similar but stronger anomalous light scattering 
properties. 7 Conductors with positive permittivity do not sup- 
port CSPs either, but they can also exhibit light transmission 
anomalies. 8 ' 9 ' 10 (Extraordinary transmission through gratings 
can even occur for acoustic waves, 11 which is completely ir- 
relevant to SPs.) Because of these contradictions, the origin 
of anomalous light scatterin g fr om metallic nanostructures is 
still being argued (e.g. Refs. [l2lU4h . 

Using modern computing techniques one may numerically 
solve Maxwell's equations for various complicated structures, 
but such computations have been largely focused on the EM 
fields. Surprisingly, the detailed mechanisms of free electron 
oscillation have been almost completely ignored in the liter- 
ature although they are known to play the essential role in 
the SP picture. Very recently, we have briefly reported our 



computations of surface charge density waves (SCDWs) and 
the role they play in the process of enhanced light transmis- 
sion through slit and hole arrays. 9 In this paper, we give a 
detailed and comprehensive illustration of the basic mech- 
anism about light emission and interference from incident- 
wave-driven free electron oscillations, demonstrate that it is 
involved in light scattering from all periodic and nonperiodic 
conducting structures (including perfect conductors), and thus 
establish a simple, concrete and universal spoof SP picture. 
This picture may provide solid guidelines for designing nano- 
optical devices by directing people to concentrate on the ge- 
ometrical parameters of conducting nanostructures so as to 
control the locations, strengths, and interference of the charge 
oscillation-induced light sources. 

n. CHARGE OSCILLATION-INDUCED LIGHT 
EMISSION AND INTERFERENCE 

To illustrate the main picture, we start from the well-known 
principle of Thomson scattering of x-rays by electrons^ in 
which the incident x rays (EM waves with wavelengths ~ 0.1 
nm) force the electrons in atoms (not necessarily free elec- 
trons) to oscillate with the same frequency. According to the 
fact that accelerating charges radiate, the oscillating electrons 
then emit new wavelets, which form the scattered waves. In 
principle, this effect also exists in the long wavelength range 
(say A > 0.1 /im), where electrons still oscillate with the inci- 
dent wave (giving rise to oscillating polarization of the atoms). 
However, since now A is much larger than the atoms (~ 0.1 
nm), the net charge density averaged on the wavelength scale 
is zero in the bulk. Net polarization-induced charges do exist 
on surfaces (or interfaces), but for non-conducting materials 
where electrons are bound to atoms and cannot move freely, 
the formation of net oscillating charges is very small even on 
rough surfaces. 

A metal has free electrons, which move/oscillate easily on 
the surface in response to external EM waves and thus may 
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FIG. 1: CSP on a planar metal surface, (a) Excitation of the CSP 
by Otto geometry. The SCDW is the result of electron oscillations 
(indicated by the dashed arrows) while the positive charges are fixed, 
(b) Symmetric sub-wavelets (Ei and E2) from a period of the charge 
wave (outlined by the dashed lines) along 6 = 0. Pi and P2 repre- 
sent two oscillating dipoles with opposite directions, caused by the 
electron oscillation, (c) Near fields of the CSP (while far fields tend 
to zero along any direction). 
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FIG. 2: Incident-wave-driven electron movement on (a) a flat con- 
ducting surface with no net charge, (b) a conducting surface with a 
single slit (net charges and their oscillation at the slit corners giving 
rise to a radiative light source), and (c) a conducting surface with a 
subwavelength slit array, (d) Periodic net charge pattern on the up- 
per surface of (c). All the dipoles and electric fields have a common 
oscillating factor e ! "* (t the time). 



emit new wavelets. But first note that a CSP corresponds to a 
surface-bound mode on the metal. If the oscillating charges 
emit light, how can the CSP be non-radiative? To clarify 
this ambiguity, let us see the Otto geometry in Fig. Q2a) as 
an example. 6 At a specific incident angle sp [greater than the 
critical angle arcsin(l/n p ) of the prism-vacuum interface], 
the incident wave can excite a CSP, which is a sinusoidal 
SCDW on the metal surface with a wavevector 

fc sp = K(l + l/e c )- 1/2 , (1) 

where e c is the permittivity of the metal and K = 2tt/\{\ 
the incident wavelength in vacuum). Here 9 sp must satisfy 
n p K sin 6* sp = Re(fc sp ), where n p is the refractive index of 
the prism. Under this condition, the incident energy is largely 
transferred to the CSP, giving rise to a reflection dip, as can be 
proved by Fresnel theory £ 

Note that CSPs can be activated only on metals with 
Re(e c ) < [and meanwhile Im(e c ) being small]. 6 The rea- 
son is that under this condition, the spatial period of the CSP 
satisfies 

d = 27r/Re(fc sp ) < A (2) 

based on Eq. (UJ. Therefore, the CSP is a subwavelength 
charge pattern compared with the incident wavelength A. 
Consider each period of the CSP in Fig. [Tfa) as a scatter unit 
that emits new wavelets. Along any arbitrary direction 9 7^ 0, 
the wavelets emitted from two adjacent units have a path dif- 
ference 

AS = dsinO < d < A, (3) 



i.e., the phase difference is less than 27r. This means that the 
oblique wavelets can never be in phase. Thus, they tend to 
cancel each other out in the far fields. The wavelets along the 
vertical direction 6 = 0, however, are in phase (AS = 0), 
but viewed from a single period of the sinusoidal SCDW 
[Fig. |TJb)], each wavelet consists of two sub-wavelets with 
opposite electric fields Ei and E2 that are also cancelled in 
the far fields. Therefore, all the emitted wavelets cannot es- 
cape the surface along any direction, so they form evanescent 
waves near the surface [Fig. \V[c)]. This gives a simple pic- 
ture why a CSP corresponds to a surface-bound mode. The 
CSP can thus propagate outside the prism-covered region in 
Fig. [Ha) without radiation loss, and the propagation distance 
solely depends on the absorption of the metal (Ohmic loss). 
Here it is obvious that media with \e c \ — > 00 (for perfect con- 
ductors) or Re(e c ) > do not support CSPs as the wavevector 
in Eq. (UJ cannot satisfy Eq. (|2]i under these conditions. 

Now we consider in Fig. |2fa) a plane wave incident on a 
conducting surface without the prism. For normal incidence, 
the incident electric field Ej n drives free electrons on the sur- 
face to move homogenously. So there is no net charge, and 
the reflection obeys Fresnel equations. 16 In Fig. |2jb), a slit 
(or hole) is added. Apparently, the electron movement now 
can be impeded near the slit corner. Here some electrons may 
move continuously to the vertical slit wall, but such move- 
ment corresponds to a 90° deflection (acceleration) while the 
incident wave does not directly provide the necessary large 
driving force. So it is reasonable to assume that most of the 
moving electrons are stopped near one corner, while positive 
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charges appear at the opposite corner because some electrons 
have moved away. This leads to the formation of an elec- 
tric dipole P a at the slit opening. P a oscillates with the in- 
cident wave with a time factor e luJt (u the angular frequency 
of the incident wave), thus acting as a new light source emit- 
ting wavelets. Such a process is in fact a Thomson scattering 
process in the optical frequency range. 

Next we apply this process to the one-dimensional (ID) pe- 
riodic slit array in Fig. |2jc). For simplicity, we assume that 
the grating is semi-infinite so that there is no feedback from 
below. Similar to Fig. |2|b), now each slit becomes a light 
source, but along any oblique direction 6^0, the wavelets 
emitted from two adjacent sources have a path difference 
AS = dsmd, where d is the period of the slit array. For an 
incident wavelength A > d, Eq. <(3j is satisfied again. Then the 
oblique wavelets are cancelled out in the far fields (destructive 
interference, similar to the absence of x-ray diffraction at non- 
Bragg angles), i.e., they also form evanescent waves near the 
surface. (This principle can also explain the fact that no light 
diffraction occurs from single crystal lattice, where the lattice 
constants are much smaller than the wavelength although the 
electrons still oscillate with the incident wave.) 

The charge pattern in Fig.|2|c) is similar to the CSP picture 
in Fig. Q] i.e., they are both subwavelength charge patterns 
(d < A). However, there are two distinct differences. First, 
the CSP is a propagating wave with a specific wavevector de- 
termined by the metal's permittivity e c in Eq. ((T), while the 
charge pattern in Fig. |2jc) is a standing wave (but not sinu- 
soidal) with the period always equal to the grating period d. 
So the former is an intrinsic property of the metal (depending 
on e c ) while the latter is a geometrical effect that can occur 
for any incident wavelength A > d and for any conducting 
materials containing free electrons [including perfect conduc- 
tors and conductors with Re(e c ) > 0]. Second, as mentioned 
above, a CSP is a complete surface-bound mode. In contrast, 
the oscillating charge pattern in Fig. |2J C ) is radiative along 
6 = (AS = 0). This can be seen from Fig. Etd), where 
we have discarded the oblique evanescent wavelets and added 
the dipoles P r that are ignored in Fig. [2c). In addition to the 
wavelet Ep a emitted from P a , P r also emits a wavelet E r 
along 6 = with a phase that is usually very close to that of 
the Fresnel reflected wave. So here we let E r include Fresnel 
reflection for convenience in discussions. Then the wavelet 
emitted from a period consists of two sub-wavelets Ep a and 
E r along 6 = with opposite directions (phases). But un- 
like Fig.QJb), Ep a and E r generally have different strengths, 
so they cannot completely offset each other. This leads to a 
propagating backward wave. Therefore, the charge pattern il- 
lustrated in Figs. |2c) and |2jd) is not a CSP, but one might 
call it a spoof SP due to its similarities to the true CSP in 
Fig.Ula)^ 

From Fig. [2] it is not difficult to obtain a general picture 
about light scattering from structured (or rough) conducting 
surfaces (either periodic or nonperiodic). When light is in- 
cident on a non-planar conducting surface, it drives the free 
electrons to move, but the movement can be impeded by the 
rough parts (particularly sharp edges) of the surface to form 
inhomogeneous oscillating charges, which become new light 



sources to emit wavelets. It is the interference between these 
wavelets that may give rise to anomalous reflection or scatter- 
ing. In the following, we will numerically prove this mech- 
anism in the simple and well studied case of periodic ID 
gratings using the rigorous coupled-wave analysis (RCWA) 
technique. 



III. RCWA OF ID LATTICE 

For monochromatic waves in a nonmagnetic medium (per- 
meability p, = 1), the electric and magnetic fields are coupled 
by Maxwell's equations (in c.g.s. units) 

V x E = -iKH, (4) 
VxH = iKsE, (5) 

where K = 2ir/\ and e is the effective permittivity. The ef- 
fective permittivity of a conductor can be expressed as e c = 
e' c — iAira/ui, where e' c is the regular permittivity and a is 
the conductivity. 16 For perfect conductors, a — » oo so that 
Im(e c ) — » — oo (which can also be derived from the Drude 
model of electrical conduction). The divergence of Eq. (O 
gives V • (eE) = 0, or 

V-E=-[(Ve)-E]/e = 47r/o, (6) 

where p(r) is the bulk charge density (including both free and 
polarization-induced charges). In a modulated medium with 
varying e(r), Ve ^ 0, which generally leads to inhomoge- 
neous charge densities p(r) according to Eq. ©. Mathemat- 
ically, e is discontinuous across a sharp interface (i.e., Ve — > 
oo), so one has to use the surface charge density Airp s = 5E Z 
to describe the charge distribution on the interface, where 
8E Z is the jump of the perpendicular electric field component 
across the interface. As an exception, TE-polarization in a ID 
structure satisfies p = since (Ve) • E = 0, so here we ignore 
it* 

From Eqs. © and one can obtain a second-order dif- 
ferential equation V x (e _1 V x H) = K 2 H., of which the 
Fourier transformation form for the ID lattice in Fig.[3]is 

K H ?n = S ' Cm — n^-m X k T1 X H n (7) 

n 

based on the Fourier expansions e _1 (x) = Cm e ~ mX and 
H(r) = X)H m e~ lkm ' r , where G m = 2nm/d (to being in- 
tegers), Cm = d' 1 f s~ 1 (x)e tGmX dx, k m = k + G m x, 
and ko is the forward wavevector. Eq. ^ can be numerically 
solved by RCWA. Here we briefly mention its main principles. 

In Fig. [3] the incident wave is Hi n exp(— «K; n ■ r) with 
Kj n = K (sin #x + cosdz). The forward wavevector ko can 
be written as ko = k$ x ~x. + qz with fcoz = Ksind, where 
q is to be determined by the eigenequation. Then the inter- 
nal diffracted wave vectors have the form k m = fc mx x + qz 
with k mx = ko x + G m . Each diffraction order to cor- 
responds to two diffracted waves exp(— iK^ ■ r) and 
exp(— iK.f n ■ r) above and below the grating, respectively. 
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FIG. 3: Geometry of light diffraction from the ID periodic grating. 
The two slit walls are located at x — ± W/2 (plus any multiple of d). 
The vertical components of the internal wavevectors k m are gener- 
ally complex vectors, x, y, z are unit vectors along the x, y, z axes, 
respectively. 



Based on the conservation of the tangential wavevector com- 
ponents, we have 



K T =K R 

mx mx 



K T = 

mz 



-K R = 

mz 



{K 2 -k 2 mx ) l l 2 for \k mx \<K, 
—i(km X — K 2 ) 1 / 2 otherwise. 



(8) 



Here note that when \k mx \ > K, the corresponding external 
waves become evanescent waves H R T e~ lkm * x exp[— {k 2 nx — 
if 2 ) 1 / 2 ^!] along ±z. In particular, for normal incidence 
(/cox = 0) and A > d, all the external waves except for 
to = are evanescent, H R ' T e~ 2 ™ lmx / d exp[—2n(m 2 /d 2 — 
1 / A 2 ) | -s |] . This is the mathematical description of the 
evanescent EM waves of the spoof SP described in Fig. [2] 

For TM polarization, all the magnetic fields are parallel 
to y. If we retain 2M + 1 diffraction orders (0, ±1, • • • , 
±M), Eq. © can be written as a (2M + 1) x (2M + 1) ma- 
trix eigenequation. From this eigenequation and the boundary 
conditions (continuity of the tangential electric and magnetic 
fields) at the two surfaces z s — 0, r, one obtains 4M + 2 sets 
of eigenvalues q. t and eigenmodes {H 7 l n } (m — —A/, — M+l, 
• • • , M and i = 1, 2, • • • , AM + 2) inside the grating and two 
sets of external fields {H*} and {ffjj (see Refs.0 and[20| 
for details). Then the zero-order reflectivity and transmissiv- 
ity are R = \H^/H in \ 2 and T = \H^/H in \ 2 , respectively. 
Meanwhile, the electric fields E r J n , E^j and (parallel to the 
xz plane) are also obtained from (the Fourier transformation 
of) Eq. (|5]l. Then the bulk charge density inside the grating 
(0 < z < t) can be calculated from 
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(k^-E^)exp(-zk^-r), (9) 



Aitp s (x) = -E h 



q z z. The surface charge density is 
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on the lower surface z = r, where 4>i = e * 9lT and = 
E^ nz exp(— iK^r). [For large r, one may need to make 
the substitutions H m e- iq ' r -> ff^ and # TO -> H^e iq ' r for 
Im(g 4 ) > to avoid numerical overflow in computing e~ iq,T , 
where H^ n is the corresponding internal wave amplitude at the 
lower surface.] For a semi-infinite grating (r — » oo), only half 
of the eigenmodes with Im(g l ) < are valid, so we only need 
to use the boundary conditions at the upper surface to com- 
pute the reflectivity and charge densities. Overall, RCWA is a 
first-principle method with the computation precision only de- 
pending on the number of diffraction orders (2M+ 1) retained, 
but note that it usually converges much slower in calculations 
of charges and near fields than in calculations of (far-field) 
reflectivity and transmissivity. 

In calculating the bulk charge density p(x, z) using Eq. (O, 
we found that when a large number of diffraction orders are 
retained, p(x, z) approaches a delta function across the walls, 
which means that "bulk" charges only exist on the slit walls, 
i.e., they are also surface charges. 9 Mathematically, we let 
p(x = ±W/2, z) represent the surface charge densities on the 
slit walls (in arbitrary units). With sufficient orders retained, 
this approximation does not affect the shapes and phases of 
the real surface charge density curves on the slit walls. 



IV. SEMI-INFINITE GRATINGS 

In the above RCWA descriptions of light scattering from a 
ID (or 2D) lattice, the eigenmodes form pairs, each pair con- 
sisting of two eigenmodes with opposite (complex) vertical 
wavevectors q and —q, respectively. As will be demonstrated 
later, one mode propagating along — z corresponds to reflec- 
tion from the bottom surface for finite r. This mode can res- 
onate with the forward-propagating one (along +z). In order 
to verify the picture in Figs.|2|c) and|2jd) without the com- 
plication of the resonance, we first consider a semi-infinite 
grating (r — » oo) where the backward eigenmodes do not ex- 
ist. 

Figure |4ja) shows the reflectivity curve calculated with 
RCWA from a semi-infinite gold grating (practically r > 200 
/xm) under normal incidence (with the frequency-dependent 
permittivity data of gold taken from Ref. 21). As a reference, 
the dashed line is the (1 — W/d)Rf curve with Rf being the 
Fresnel reflectivity from a flat gold surface (Rf — 0.98 in the 
wavelength range 0.8 — 10 pm) and W the slit width. Com- 
pared with this reference curve, the anomalous reflection phe- 
nomenon from the grating is obvious. Generally the reflectiv- 
ity Rq is less than (1 — W/d)Rf except that near the Wood's 
anomalies A ~ d/|m| (m ^ being integers), Rq is close to 
unity. 

FigureUb) shows the charge density function p s (x) on the 
upper surface (z — 0) for an arbitrary wavelength in the A > d 
range. This function correctly shows that the incident wave 
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FIG. 4: (Color online) Reflectivity (a) and surface charge densities 
(b) of a semi-infinite gold grating with period d — 2 /im and slit 
width W = 0.2 ^m. Normal incidence. 



indeed causes significant inhomogeneous charges on the grat- 
ing surface with the charges highly accumulating near the slit 
corners, which is excellently consistent with the charge dis- 
tribution pattern predicted in Fig. |2jd). At a time when E; n 
is toward +x at z = 0, we have predicted in Fig. |2d) that 
the phase of the charge pattern is constant, equal to — 7r (neg- 
ative charges), on the left half surface W/2 < x < d/2, while 
for d/2 < x < d — W/2, the phase is (positive changes). 
FigureHJb) shows that this prediction is largely correct except 
that the calculated phases are slightly displaced from the pre- 
dicted phases and it by Aip ~ 0.17T in most regions on the 
surface. The phases near the slit corners are closer to the pre- 
dicted values. Our calculations show that the charge patterns 
are nearly the same for any wavelength A > l.ld with no res- 
onance, and the phase shift Aip decreases with increasing A, 
i.e., Aip — > for A 3> d. Therefore, the calculations indeed 
confirms the picture of charge accumulation and oscillation on 
the subwavelength lattice in Figs.|2|c) and[2fd). Apparently, 
the period of the charge pattern in Fig.Ub) is strictly equal to 
the lattice constant d and is irrelevant to the dispersion prop- 
erty of CSPs in Eq. (Q~|i. By performing RCWA calculations on 
gratings made of conductors with Re(e c ) > or perfect con- 
ductors with Im(e c ) — > — oo, we found that the main features 
of the charge patterns remain unchanged. 

To understand the anomalous reflection for A > d in 
Fig.HJa), we may simply consider that, as we mentioned be- 
fore, the E,. wavelet in Fig. |2|d) consists of two contribu- 
tions, E r = —(Ef + E' r )5t, with Ef corresponding to regu- 
lar Fresnel reflection (Ef = —y/RjE^ and E' r correspond- 
ing to the emission of light from dipole P r along the back- 
ward direction 9 = 0. Without charge accumulation, we have 
Ep a = E' r = and the reflection should obey the Fresnel 



theory Rq = (1 — W/d)Rf. When charges appear at the 
slit corners, it can be verified by RCWA that P a is strength- 
ened much faster than P r (for W < d/2). Then the effective 
strength of Ep a becomes larger than E' r . Consequently, Ep a 
completely cancels E' r , and also partially offsets Ef. Thus, 
the net effect is that the overall reflectivity Rq is smaller than 
(1 - W/d)R } . 

When A is reduced to be less than (or close to) the grating 
period d, some of the external wavevectors K^ z in Eqs. ([8]) 
become (or tend to be) real, and the corresponding diffracted 
waves become non-evanescent. Then the diffraction effect ap- 
pear, which can significantly changes the reflectivity (partic- 
ularly at the Wood's anomalies A ~ d/\m\). The details in 
the diffraction range A < l.ld are discussed in the Appendix 
since they no longer belong to subwavelength optics. But, it 
is worth emphasizing again here that the diffraction effect is 
absent for the entire long wavelength range A > l.ld, where 
the subwavelength charge patterns are always nearly the same 
as that in Fig. @Ib). 



V. FINITE-THICKNESS GRATINGS 

In Fig. |2] we have indicated that the dipole P a also emits a 
wavelet in the slit toward +z [see E a in Fig. [5ja), which may 
also include a portion of the incident wave directly transmitted 
into the slit]. Due to the waveguide constraint, E a tends to 
be a plane wave inside the slit, i.e., E a ~ E a exp(~ik z z)x, 
where k z ~ 2ir/\. Similarly, this wave drives electrons on the 
slit walls to oscillate, resulting in two SCDWs p a exp(—ik z z) 
and — p a exp(—ik z z) (with p a oc E a ) on the two opposite 
walls, respectively. The SCDWs and the E Q wave propagate 
along +z and attenuate gradually due to the absorption of the 
conductor. If the grating is extremely thick, these waves can 
be completely absorbed before reaching the bottom surface, 
which corresponds to the semi-infinite case. 

If the grating is thin enough, the SCDWs on the walls can 
reach the exit surface without significant absorption. Then 
in a similar way, the moving charges can be impeded at the 
lower slit corners, leading to another large oscillating dipole 
Pi,, as shown in Fig. [5ja),- P h can give a strong feedback to 
the upper surface by emitting a wavelet Ef, ~ E exp(z£; z z)x 
propagating upward. Ef, also corresponds to two SCDWs, 
±p exp(ik z z), on the two walls, which are the back-bounced 
SCDWs of the ±p a exp(— ik z z) waves by the bottom corners. 
If E;, is in phase with E a at z — 0, it enhances P a . The en- 
hanced P Q subsequently strengthens Ep a , E a , P;,, Eh, and so 
on. Then a Fabry-Perot-like resonant state is formed, with E Q 
and Eb forming a standing wave E a + Eh oc cos{k z z)e lut in 
the slit. Under this condition, E r is largely offset by Ep Q in 
the far fields, leading to minimized backward reflection. Fig- 
ure|5jb) shows the zero-order reflectivity curve of a gold grat- 
ing with thickness t = 4 p,m. Compared with Fig. H|a), one 
can see that Figure [2b) indeed shows a number of reflection 
dips corresponding to Fabry-Perot resonance. 

At the exit surface z = t [Fig. [5ja)], dipoles P^ and P t 
also emit wavelets towards the outside of the slit. For A > d, 
only the wavelets Ep;, and E t can propagate along +z (while 
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FIG. 5: (a) Process of charge oscillation-induced light emission, res- 
onance and tranmssion through a conducting grating. Reflectivity 
(b) and transmissivity (c) of a gold grating with parameters d = 2, 
W = 0.2, and r = 4 pm. 



the oblique wavelets again form evanescent waves). Unlike 
the case above the upper surface where E r contains specular 
reflection, here wavelet Et is purely emitted from dipole Pt. 
For W < d/2, the strength of (Epb) is much stronger than 
that of P t (Et), so the transmitted wave is dominated by Ep;,. 
Consequently, the energy of the transmitted wave is highly 
localized near the exit opening. For long wavelengths A ^> 
d, such a "near-field focusing" effect can achieve a focusing 
width W far smaller than A, which has potential applications 
in nano-focusing/beaming, lithography, etc. In the far-field 
region, however, this effect disappears as the transmitted beam 
becomes a plane waveM 

At resonant wavelengths, since the strength of wavelet 
Ept is maximized, the zero-order transmissivity T is also 
maximized, 9 as can be seen in Fig.|5jc), where each reflection 
dip exactly corresponds to a transmission peak (also see sim- 
ilar results from finite-difference time-domain calculations in 
Ref-Q. 

If the waves E a and Ef, are ideal plane waves with wavevec- 
tor k z = 2ir/\, Fabry-Perot resonance should occur at Ajv = 
2t/N (N the resonance order) except that resonant peaks 
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FIG. 6: (Color online) Charge patterns on the thin gold grating, (a) 
Charge densities on the two surfaces z s = and r for the resonant 
wavelength A3 = 3.19 /J,m of peak TV = 3 in Fig. [5] (b) Nearly 
standing charge wave on the slit wall x = —W/2 at A3. Note that 
p v (W/2, z) = — p v (— W/2, z). (c) Surface charge densities at a 
non-resonant wavelength A = 3.76 /im (corresponding to the valley 
between peaks N — 2 and N = 3 in Fig.|5]l. (d) Charge wave on the 
slit wall x = — W/2 for A = 3.76 fim, where z — is no longer an 
anti-node of the standing wave. 



with Ajy < d are suppressed by the diffraction effect. 9 How- 
ever, the actual resonance wavelength is always redshifted, 
Xn = 2t/N + An, where the redshift An may vary (slowly) 
with d, W, t, and e c . One reason for the redshift is that the 
standing wave is distorted near the two ends of the slit [see 
Fig. |6jb)]. One may refer to Ref. IH for discussions of other 
possible mechanisms. Here note that due to the redshift, the 
spatial period of the (sinusoidal) SCDWs on the wall is less 
than the incident wavelength by approximately Ajv, i.e., they 
are also subwavelength charge waves. Based on An ~ 2t/N, 
thick (and highly conducting) gratings (r 3> d) have many 
resonance wavelengths in the non-diffraction range A > d, as 
experimentally demonstrated in Ref. |24j. 

In Figs. Oa) and Ifjb), the computed charge density dis- 
tributions at resonant wavelength A3 are well consistent with 
the picture of Fig. [5J a )- in Fig- 13 a X the two surface charge 
patterns \p s (x, z = 0, r)| are very similar to that in Fig. |4{b), 
which confirms the existence of the large dipoles P a and Pf, in 
Fig. |3 a). Note that the charge densities in Figs.|6ja) andHfb) 
are in the same (arbitrary) unit. Therefore, the charge densi- 
ties near the slit corners are much higher in Fig. |6|a) than in 
Fig. |4fb) due to the Fabry-Perot resonance/enhancement. As 
also shown in Fig.[6{a), in thin gratings where the attenuation 
of the charge density waves on the slit walls is negligible, the 
two SCDWs p s (x, z = 0,t) are almost identical except that 
for odd resonant orders TV, they have a phase difference ir. For 
relatively thicker gratings, the strength of p s (x, z = t) drops 
with increasing r. For r > 200 p,m, p s (x,z — r) almost 
disappear while p s (x, z = 0) tends to be the same as that in 
Fig.Hb). 
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Figure |6jb) correctly reveals that on the slit walls, the 
charge density waves p(x = ±W/2, z) with approximately 
stepped phases are nearly standing waves. Here the p pro- 
file also shows high accumulation of charges at the slit cor- 
ners (x = ±W/2, z = 0, t) that are (always) in phase with 
p s {x = ±W/2, z = 0,t)£ 

In Fig. [3J a), if Eh not in phase with E Q (and E; n ) at z = 0, 
it suppresses the strengths of P Q and Ep a and influences their 
phases. Consequently the strengths of the charge waves on 
the slit walls are also reduced, leading to a weaker dipole Pf, 
and weak transmissivity. This mechanism is clearly shown in 
Figs.|6jc) and[6fd) at a non-resonant wavelength. Compared 
with Figs.|6ja) and|6fb), the charge densities at the slit corners 
(x = ±W/2, z = 0, r) all drop significantly for both p and 
p s , particularly at the upper corners. Meanwhile, the phases 
of the charge waves are also altered so that no resonance is 
formed. As stated above, without surface charges, the reflec- 
tivity from the upper surface should be the Fresnel reflectivity 
Rq = (1 — W/d)Rf. Here one can see from Figs.|6jc) and 
[6fd) that at non-resonant wavelengths, the strengths of P a and 
P r at the upper surface are very small, and then the reflectiv- 
ity Ro in Fig.|5jb) is indeed very close to (1 — W/d)Rf ~ 0.9 
in most of the non-resonant wavelength range. For the same 
wavelength, the non-resonant reflectivity in Fig.[5jb) is much 
stronger than that in Fig.@fa) where charge oscillation is heav- 
ily involved. This further proves the essential role charge os- 
cillation plays in extraordinary light scattering from metama- 
terials. 

As mentioned above, perfect conductors with Im(e c ) — > 
— oo and conductors with Re(e c ) > do not support CPSs. 
However, we have demonstrated in Ref. 9 that ID gratings 
with Re(e c ) > may still show similar extraordinary light 
transmission although the transmissivity is relatively lower. 8 
Here we use RCWA to simulate the transmission through 
a nearly perfectly conducting grating with a large constant 
imaginary permittivity e c = — ilOO 000. Based on this value, 
the wavevector fc sp in Eq. (Q~|i is almost accurately equal to K, 
so CSPs should not exist. However, our calculations show that 
all the major properties in this case are almost identical to that 
of regular metallic gratings. For example, Figure Q a) shows 
the transmissivity curve calculated with the same geometrical 
parameters in Fig. |3c), while Figure |7Jb) shows the charge 
densities on the grating surfaces for the third-order resonance 
peak. Compared with the reflectivity curve in Fig. [51c) and the 
charge density distribution in Fig. |6ja), Figure Q apparently 
indicates that the light scattering mechanisms for the perfect- 
conductor case are the same, and thus are irrelevant to CSPs. 
In fact, the resonant transmissivity peaks and the charge den- 
sities in Fig. [7] are averagely higher than those for the gold 
grating, indicating that high conductivity can (significantly) 
enhance the extraordinary scattering effects. 

As demonstrated in refs.l9landl25l extraordinary transmis- 
sion or scattering through 2D hole arrays involve the same 
mechanisms of light emission and interference except that the 
tunneling of the SCDWs through the holes is different. The 
details of oblique incidence geometry will be presented else- 
where, but the charge oscillation principle is similar. 
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FIG. 7: (Color online) Light scattering from a nearly perfectly con- 
ducting grating. e c = — ilOO 000 (compared with e c = —2834 — 
il333 for gold at A = 9.9 pm). d = 2, W = 0.2, and r = 4 
pm. Normal incidence, (a) Transmissivity curve. Compared with 
Fig-H c X h ere the redshifts (Ajv) of the resonant wavelengths are 
smaller (so peak N — 5 below A = d is truncated) and the resonant 
transmissivities are larger. 



VI. NONPERIODIC STRUCTURES 

From the above demonstrations, it becomes obvious that 
charge oscillation-induced light emission and interference 
are a fundamental and universal mechanism underlying var- 
ious extraordinary light scattering processes from conducting 
structures although these processes may also involve other 
mechanisms simultaneously (e.g. cavity resonance). The 
only basic requirement for this mechamism to work is that 
the structure have free electrons. So this mechanism applies 
to structures of metals, perfect conductors, conductors with 
Re(e c ) > 0, semiconductors, etc, but high conductivity can 
significantly enhance the anomalous scattering effects. 

This mechanism also applies to nonperiodic structures.— 
From Fig. |2|b) one can see that an isolated single slit also 
acts as a light source. If the conducting plate has a finite 
thickness, the exit opening of the slit at the lower surface be- 
comes another strong light source at a Fabry-Perot resonant 
wavelength, emitting a transmitted beam below the plate. 23 
Compared with the periodic slit array in Fig. |5Ja), the waves 
emitted from a single slit have no interference and thus are 
completely divergent. If the slit is surrounded by periodic 
grooves on the entrance surface, as shown in Fig. [8la), each 
groove now acts as a light source. Under the conditions that 
the groove period is less than the incident wavelength and 
that the Fabry-Perot-like resonance can be achieved simulta- 
neously in both the grooves and the slit, the EM fields above 
the upper surface become similar to those in Fig. |5ja) with 
the oblique waves forming evanescent modes. Most impor- 
tantly, the enhanced backward wavelets from the light sources 
[Ep a in Fig. |5ja)] significantly reduce the Fresnel reflection 
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(E r ). Accordingly, the backward reflection is reduced while 
the transmission through the slit can be greatly enhanced (by 
up to 2 orders in Ref. |27|) . But, the transmitted beam below 
the plate is still divergent. 

Now if similar grooves are made on the exit surface, the 
wavelets emitted from the exit opening of the slit also drives 
free electrons to form oscillating dipoles at the openings. 
Thus, they also become light sources. Similarly, the oblique 
wavelets emitted from the slit opening and the grooves on the 
exit surface tend to form evanescent waves, giving rise to a 
narrow and directed transmitted beam below the slit opening. 
However, the dipoles on the upper entrance surface are formed 
and driven mainly by the wide incident wave while those on 
the exit surface result only from the wavelets emitted from the 
dipole of the single slit. Therefore, the strengths of the light 
sources on the exit surface decrease quickly with increasing 
distances of the grooves from the slit. Due to this reason, 
the transmitted beam cannot be completely collimated since 
the oblique wavelets cannot be completely suppressed. Mean- 
while, the grooves on the exit surface have little influence on 
the overall transmission efficiency^ 

Obviously, the picture illustrated in Fig. [HJa) can also ex- 
plain enhanced light transmission and directed nanobeaming 
through a single aperture surrounded by circular grooves in 
the 2D case except that the cavity resonance mechanisms in 
the aperture and in the circular grooves may be different and 
the directions and distributions of the light sources near the 
groove edges are more complicated^ 

As another nonperoidic structure example, it is known that 
when one end of a conducting nanowire is illuminated by a 
narrow-wavefront beam (with the electric field E[ n being po- 
larized along the wire), the other end that is not illuminated 
can emit light, as schematically illustrated in Fig. [8jb). The 
common explanation of this phenomenon is that light is trans- 
ferred by CSPs on the wire surface. 28 However, it is found that 
this phenomenon is more pronounced in the low-frequency 
(e.g. terahertz) range, where most metals become nearly CSP- 
free perfect conductors. In particular, Wang and Mittleman 7 
have experimentally demonstrated that in the terahertz range 
(A ~ 1 mm), the wave modes on metallic nanowires have a 
dispersion trend that is opposite to that of CSPs. 

In fact, according to our charge oscillation picture, the ba- 
sic mechanism underlying light transfer on nanowires is very 
simple. In Fig. |8jb), the incident wave drives free electrons 
near the left input end to oscillate. The agitated electrons then 
propagate outside the illumination area towards the right side 
as a SCDW. At the other end the propagating charge wave 
is discontinued, giving rise to a strong charge accumulation 
there. The oscillation of these charges them emit new light 
near the exit end. Meanwhile, the charge wave is bounced 
back. When the bounced charge wave is in phase with the 
forward wave (and the incident wave) at the input end, Fabry- 
Perot resonance occurs. This is very similar to the charge 
movement on the slit walls in Fig. EJa). In general, the res- 
onant wavelength here also has a redshift. Accordingly, the 
standing charge wave on the wire surface is a subwavelength 
wave (d < A), and based on Fig. [T] it generates little radia- 
tion loss when propagating on the wire. The Fabry-Perot res- 
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FIG. 8: (a) Schematic representation of enhanced light transmis- 
sion and directed nanobeaming through a single slit surrounded by 
grooves, (b) Transfer of light on a conducting nanowire by subwave- 
length SCDWs. 



onance and the subwavelength charge patterns (proportional 
to the strengths of the near fields) indeed have been demon- 
strated both experimentally and theoretically . 7 ' 28 ' 29 . Interest- 
ingly, our calculations show that a charge wave propagating 
in the unilluminated region of a flat/straight conducting sur- 
face (including the slit wall and the straight wire) is always 
a subwavelength SCDW, which indicates that the general SP 
picture elaborated in the literature, light — > subwavelength 
SCDWs — > light, is indeed correct except that the SCDWs are 
spoof SPs and do not necessarily have the dispersion property 
ofEq. ©. 

According to this picture, the conductivity of the nanowire 
is the dominant factor determining the propagating distance 
of the charge waves and the efficiency of light transfer. This 
explains the remarkably high transfer efficiency in the long- 
wavelength range where most metals are highly conducting. 
To further confine the near fields so as to reduce the radia- 
tion loss (caused by possible deviations of the actual charge 
waves from ideal subwavelength standing waves), one may 
activate charge waves in grooves and guide them to propagate 
inside the grooves. In these cases, the charge waves are chan- 
nel spoof SPs 30 that may have longer propagating distances. 

From Fig. Ob) it is obvious that to achieve high transfer ef- 
ficiency, the diameter (vertical dimension) of the wire should 
be much smaller than the incident wavelength so that the ag- 
itated charge waves on the top and bottom of the wire have 
nearly the same phase. Otherwise, the charge waves with dif- 
ferent phases will quickly mix together and thus offset each 
other outside the illuminated region, leading to a short propa- 
gation distance. This explains why light transfer is remarkable 
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on nanowires. With respect to this effect, it is expected that a 
thin conducting slab would be more efficient since it can en- 
hance the input coupling efficiency and reduce the electrical 
resistance without causing phase differences. 

Note that in Fig. |2c), when the moving electrons are 
stopped at the slit corners, they also have a tendency to 
be bounced back, similar to the moving electrons on the 
nanowires in Fig.[HJb) [and on the slit walls in Fig.|5ja)]. The 
difference in Fig. |2jc) is that the bounced back charges are 
suppressed by the incident electric field Ej n since the driv- 
ing force provided by Eu, is always opposite to this tendency, 
while on the nanowire of Fig. 13b), E; n is absent except at the 
input end. 



VII. SUMMARY 

By numerically calculating the SCDWs on gratings, we 
have demonstrated that the incident wave can drive free elec- 
trons to accumulate and oscillate near the slit corners to form 
new light sources. These light sources then emit new wavelets. 
For periodic subwavelength structures (d < A), the oscillat- 
ing charges form subwavelength charge patterns (i.e., spoof 
SPs) and the wavelets emitted from them destructively inter- 
ference with each other to form evanescent wave modes near 
the surfaces. Usually combined with other mechanisms (e.g. 
Fabry-Perot or cavity resonance), the spoof SPs can lead to 
anomalous light reflection, transmission, or scattering. The 
spoof SPs are mainly a geometrical effect and generally do 
not have the dispersion properties of CSPs. Note that in the 
literature, the SP-like modes on metamaterials with finite con- 
ductivity were widely assumed to be CSPs, while only those 
on perfectly conducting structures were believed to be spoof 
SPs. Here we have demonstrated that they are all spoof SPs. 
(For transmission of acoustic waves through gratings, 1 ^ the 
counterpart of charge oscillation is the mechanical vibration 
of the structured medium, particularly near the corners and 
edges, that emits acoustic wavelets.) 

We also illustrated that the same mechanism of charge 
oscillation-induced light emission and interference applies to 
all structures with free electrons (including perfect conduc- 
tors and nonperiodic structures). Thus, the spoof SP picture 
represents a basic and universal mechanism of light scattering 
from conducting nanostructures. The guideline provided by 
this mechanism is that in designing novel nano-metamaterial 
devices, there is no CSP excitation constraint, but one needs to 
precisely design the geometrical parameters of the devices so 
as to accurately control the locations of the new light sources 
(including maximizing the strength of various resonance pro- 
cesses if involved) and their interference. Meanwhile, choos- 
ing highly conducting materials is generally another require- 
ment for enhancing the anomalous scattering effects. 
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APPENDIX: CHARGE PATTERNS IN THE DIFFRACTION 
RANGE FOR ID GRATINGS 

One might intuitively think that for normal incidence with 
the electric fields Ei„ parallel to the surface in Fig. |2fc), the 
charge pattern on the surface should remain the same even 
if A < d, as Ejn provides the same driving force. How- 
ever, the free electrons on the grating surface are driven by 
both the incident wave and the wavelets emitted from the 
charges/dipoles. The latter may significantly influence the 
charge distribution (to achieve self-consistency of the system) 
when A < d. 

Figure [9] shows the charge patterns in the short- wavelength 
range (A close to or below d) for the semi-infinite grating 
of Fig. Eta). To understand these patterns, recall in Sec. [Ill] 
that the ±mth-order wave components (m > 0) above the 
grating can be written as E^ l e~ 27T " nx / d exp(—iK^ lz z) and 
E-R m g 27rima; / rf exp(— iK^ z z), respectively for normal inci- 
dence, where K^ z — — 2-7r(l/A 2 — m 2 /d 2 ) 1 ^ 2 . By symmetry, 
|E^j| = |E^ m |. So these two waves form a standing wave 
£ m sm(2Trmx/d)exp(—iK^ lz z). Accordingly, there exists 
a standing SCDW p s . m sm(2irmx/d) corresponding to this 
wave on the metal surface (but discontinued in the slit gap). 
For A > d/m, the propagating factor exp(—iK^ z z) of the 
£ m wave becomes a decaying factor exp(— /3 m |z|) along —z, 
where m = 2n(m 2 /d 2 — 1/A 2 ) 1 / 2 is the decaying coeffi- 
cient. Therefore, the £ m wave is an standing evanescent wave 
under this condition. For A > d (so that A > d/m), all the £ m 
waves except for m = are evanescent, so they all have weak 
strengths. The charge pattern in Fig.HJb) is the collective con- 
tributions from a large number of charge wave components 
p s ,m sm(2irmx/d), and this pattern (as well as the collective 
near fields of the wavelets emitted from the dipoles) coincides 
with the grating lattice. 

When A decreases toward d from above, the first-order 
£i wave tends to be strengthened as its decaying coeffi- 
cient (3i decreases. Compared with Fig. Hfb), one can see in 
Fig. |9ja) the appearance of the standing SCDW component 
p s ,i sm(2irx / d) superimposed on the overall charge pattern 
that is similar to that in Fig. Hfb). (Here the charges are in- 
fluenced by the near fields even if the £\ wave is evanescent 
in the far fields.) Meanwhile, the phase shift A(p indicated in 
Fig-Htb) increases to — 0.287T in Fig. [9^a) due to the phase of 
the complex constant p St ±. 

The wavelength A = d corresponds to the first-order 
Wood's anomaly. 9 Under this condition, the wavelets emit- 
ted from two adjacent units of the grating along the hori- 
zontal direction 8 = ir/2 have a phase difference of 2ir ac- 
cording to Fig. |2c), which corresponds to a special Bragg 
diffraction condition. Then the strength of £\ wave is max- 
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FIG. 9: (Color online) Diffraction affected charge densities on the 
surface z = of a semi-infinite gold grating, d = 2 and VK = 0.2 
/im. Dashed lines are the phases. Normal incidence. 



imized and becomes outstanding among other components, 
as can be verified by RCWA. Therefore, the charge pattern 
in Fig. |9|b) is dominated by the corresponding charge wave 
p s ,i sm(2irx/d). Here one (virtual) node (Node 1) of this 
standing SCDW is located at (x = 0, z = 0), the middle of 
the slit [while the other node (Node 2) is at (a; = d/2, z = 0)]. 
Node 1 significantly suppresses the charge densities at the two 
slit corners (because the nodes of a standing wave have zero 
amplitudes), so the strength of the dipole P a [see Figs. |2jc) 
and|2jd)] is minimized. As discussed in Figs.[6]c) and[6fd), 
when the strength of P a is very small, the reflectivity should 
approach the Fresnel reflectivity (1 — W/d)Rf. But note 



that in Fig.|9jb), the charge pattern p s> i sin(27riz;/<i) also con- 
tributes positively to the reflectivity [equivalent to the increase 
of E' r in Fig. |2d) although it has a phase difference ~ tt/2 
from the incident wave]. This is the reason why in Fig. HJa) 
the reflectivity at Wood's anomaly wavelengths is even higher 
than (1 — W/d)Rf. Note that the first-order Wood's anomaly 
in Fig.@|a) is slightly red-shifted from d = 2 to 2.01 pm. 

When A decreases to 1.9 /im in Fig.|9jc), the S\ wave is still 
a propagating mode, but its strength decreases as the diffrac- 
tion condition deviates from the Bragg condition. Meanwhile, 
the p S} 2 sin(47rx/(i) wave begins to gain strength, which mod- 
ifies the stepped phase profile. When A further decreases to- 
ward d/2 in Fig. |9jd), the p Sj 2 sm(Aivx/d) wave becomes ap- 
preciable. At the second-order Wood's anomaly A = 1.01 pm 
(also slightly red-shifted from d/2), this wave is maximized 
[Fig. Ete)]. For A < d/2 [Fig. EJT)], the third-order wave 
p Si 3 sin(67ra/d) begins to show strengths, and so on. Thin 
gratings have similar properties, but these effects are mean- 
while mixed with the resonance in the slits. 

Overall, we may call the short-wavelength range A < 1.1c? 
the diffraction range, where the diffraction effect of the ID 
lattice appear, particularly at the Wood's anomalies. In the 
diffraction process, it is interesting that the incident energy is 
largely back reflected rather than diffracted although the cor- 
responding diffracted waves become non-evanescent, as can 
be seen from Fig.Hta). This is also true in Figs. [2b) and|3c) 
for the thin grating, where the reflectivity is averagely very 
high in the entire diffraction range while the transmissivity 
(particularly at the Fabry-Perot resonant transmission peaks) 
is much lower than that in the non-diffraction range A > 1 . Id. 
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